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Abstract 



We apply the mirror principle of [L-L-Y] to reconstruct the Euler data Q = 
{Qd}denu{o} associated to a vector bundle V on CP™ and a multiplicative class b. 
This gives a direct way to compute the intersection number Kd without referring to 
any other Euler data linked to Q. Here Kd is the integral of the cohomology class b(Vd) 
of the induced bundle Vd on a stable map moduli space. A package 'EulerData_MP .m' 
in Maple V that carries out the actual computation is provided. For b the Chern 
polynomial, the computation of K\ for the bundle V = X^CP 2 , and Kd, d — 1,2,3, 
for the bundles O c p 4 (0 with 6 < I < 10 done using the code are also included. 

Key words: Atiyah-Bott localization formula, concavex bundle, Euler data, linear cr-model, S 1 x T n - 
equivariant cohomology. 
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Reconstruction of Euler Data 
0. Introduction and outline. 

Introduction. 

Ever since the ground-breaking work of [C-dlO-G-P], mirror symmetry and its meaning 
and consequenes have been investigated by several groups of people both from the math- 
ematical and from the physical point of view. (See introduction of [L-L-Y] and references 
therein for background and for a comparison of different approaches. See also [MS].) 

In this note, we apply the theory developed in the series of papers [L-L-Y] I, II, and 
III, to compute new intersection numbers on stable map moduli spaces. The theory goes 
beyond justifying mirror symmmetry that relates the usually difficult A-model computa- 
tions to the much more tractable B-model computations for Calabi-Yau manifolds. Indeed, 
the theory is a way to directly do the "A-model" computations for manifolds which are 
not necessarily Calabi-Yau. It is the goal of these notes to explain this and to provide a 
computer code that carries out the actual computations for bundles over CP™. 

To make this article more self-contained, we recall in Sec. 1 the definitions of the basic 
objects involved and the Atiyah-Bott localization formula that is used substantially in 
the theory. In Sec. 2 and Sec. 3, we focus on the case of critial bundles over CP™ and 
consider their Euler classes. In Sec. 2, we give a quick summary of facts and formula from 
[L-L-Y] I-III that are directly related to the actual computation of Euler data {Qd}d and 
the intersection numbers K^. In Sec. 3, we explain how the theory of [L-L-Y] gives rise 
to a system of linear equations that can be solved inductively. The solution of the system 
gives the Euler data {Qd}d, from which the intersection numbers can be computed. 
After these, we then discuss in Sec. 4 the modifications needed to take into account also 
non-critical bundles. There the Chern polynomial is considered. In Sec. 5, we single out 
six examples whose first few are computed this way via a Maple code. In Sec. 6, the 
Maple code EulerData_MP .m with instructions is given. Eighteen cases have been tested 
and computed. The last record of the run for each of these cases is given in SEC. 3 of the 
code for references. The code provided can be easily modified to compute other cases of 
interest. 

This article is served as a supplement to and a computational account of [L-L-Y]. As 
a result, our notations and terminologies follow [L-L-Y] very closely. Readers are referred 
to ibidem for more theoretical details. 

Outline. 

1. Essential mathematical backgrounds for physicists. 

2. Summary of related constructs in "Mirror Principle". 

3. Computation of Qd inductively. 

4. Modifications for non-critical bundles over CP". 

5. Examples. 

6. A package in Maple V for the computation of Qd and Kd- 
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1 Essential mathematical backgrounds for physicists. 



We collect in this section the most essential backgrounds for understanding these notes. 
Along the way, we also set up the notations for the notes. 

• Stable maps and their moduli. [Ko] 

Definition 1.1 [stable map]. Let X be a smooth projective variety. An n-pointed 
stable map into X consists of a connected marked curve (C,pi, • • • ,p n ) and a morphism 
/ : C — > X satisfying the following properties: 

(1) The only singularities of C are ordinary double points. 

(2) pi, ■ ■ ■ , p n are distinct ordered smooth points of C. 

(3) If Ci is a component of C that is isomorphic CP 1 and is mapped to a point 
under /, then Ci contains at least three special (i.e. nodal or marked) points. 

(4) If C has (arithmertic) genus 1 and n = 0, then / is not constant. 

Remark 1.2. Given Conditions (1) and (2) in the above definition, Conditions (3) and 
(4) are equivalent to the assertion that the data (f,C,pi, • • • , p n ) has only finitely many 
automorphisms . 

Given a class (3 G H2(X, Z), the moduli space of all stable maps (/, C,p\, ■ ■ ■ , p n ) such 
that [f(C)] = (3 into X will be denoted by M gtn (X,(3). 

• Equivariant cohomology. (See [Au].) Given a group G acting on a space X. Let BG 
be the classifying space and EG — ► BG be the universal principle G-bundle associated to 
G. The equivariant cohomology H G {X) of X associated to the G-action is defined to be 



where Xq = EG Xg X is the total space of the associated A-bundle over BG. Note that 
H^(pt) = H*(BG), and that H* G {X) is naturally a H* G {pt)-m.odu\e. 

The constant map X — > pt induces an equivariant projection ttx '■ Xq — > BG. The 
induced pushforward map ttx\ from H G {X) to H G (pt) is given by integration along the 
fiber of irx- This is also called the equivariant integral. In notation, 



Remark 1.3. In this article, the coefficient for H*{Xq) can be Q, R, or C. Usually we use 
Q or C in the discusssion. 

Example 1.4. Let T r = JJr S 1 be an r-torus. Then BT' = T7 r CP 00 and H*(BT r ) = 
Hjr(pt) = C[Ai, • • • , A r ], the polynomial ring generated by Ai, • • • , A r , where Aj is the 
first Chern class of the hyperplane line bundle 0(1) over the i th CP 00 in the product. 



H* G (X) = H*(X G ) 




2 



Let T r — > GL(N + 1,C) be a representation of an r-torus on C N+1 with weight 
((3q, /3jy). Note that each is a linear combination of the A/s. This induces a 
T r -action on CP^. With respect to this action, 

uMcp") = flf(pf)[C]/n£o(C-ft). 

where £ is the equivariant hyperplane class from a fixed lifting of the hyperplane class of 
CP^. The equivariant integral J(cp") T r ■ ^^(CP^) — > Hjr(pt) picks out the coefficient 
of C W in elements of iff (CP^). 

□ 

• The Atiyah-Bott localization formula. [A-B] Let T be an r-torus that acts on a 
manifold X with the set of fixed points a union of smooth connected submanifolds Zj. 
Then the normal bundle Nj of Zj in A" is a T-equivariant vector bundle with its equivariant 
Euler class e T {Nj) G H^(Zj). 

There are three fundamental maps between H^{Zj) and H^{X) : 

(1) the restriction homomorphism: 

i* : H* T {X) — HftZj) 

induced by the equivariant inclusion ij : Zj ^ X; 

(2) the Gysin map: 



J ! 



: H?p(Zj) — > ify(A) 



Note that ([Au]), for any aj G H^(Zj), one has 

z* o i j: (a) = aj- U e T (A i ). 

Let 1Z be the localization of H^(BT). With the notation following Example 1.4, 
7£ = C(Ai, • • • , A r ). Then ex{Nj) is an invertible element in the localization H^(Zj) <g> 1Z. 
We can now state the Atiyah-Bott localization formula [A-B]: 

Fact 1.5 [Atiyah-Bott localization formula]. The following map is an isomorphism 

H^(X) ®K © H* T {Zj) n 



a 



/is inverse is given by 

(a>j)j i — ► 51 



Combining these two, one has 



.7! 



er(A,) 
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for any a G H^(X) ® 7£. 



Example 1.6. {Continuing Example 1.4). The fixed point set on CP^ comes from the 
N + l coordinate lines of C N+1 . Denote this set by {p , ■ ■ ■ , Pn }• Let a G H£(CP N )®Tl, 
then a can be written as a polynomial /(C) with coefficients in 1Z. In terms of this, 
ij(a) = f{Pj); ex(Nj) = Y\ k _^j(Pj — /3k), where k runs in { 0, • • • , N }; and the Gysin map 
ijl is given by the cup product with Ylk^jiC ~ Pk)- Thus, from the localization formula, 
one has 

□ 

• Concavex bundles over CP n . ([L-L-Y, I].) Let T = T n+1 acts on C n+1 with weights 
Ao, • • • , A n . It induces an action on CP™. 

Definition 1.7 [concavex bundle]. Let V be a T-equivariant vector bundle over CP n . 
We call V convex (resp. concave) if the T-equivariant Euler class ct(V) is invertible and 
H\C,f*V) = (resp. H°(C,f*V) = 0) for any 0-pointed genus stable map / : C 
CP n . We call V concavex if it is a direct sum of a convex and a concave bundle. We 
denote this decomposition by V = V + © V~ . 

Definition 1.8 [splitting type]. Let V be a T-equivariant concavex bundle over CP n . 
Let l a , kb be positive integers such that for every T-invariant line C = CP 1 in CP™ we 
have a T-equivariant isomorphism 

v\ c = @ a o{i a ) © e b o{-k b ). 

Then we call ( l±, • • • ; fei, • • • ) the splitting type of V . 



2 Summary of related constructs in "Mirror Principle". 

This section follows [L-L-Y, I]. Readers should consult [L-L-Y, I] (see also [L-L-Y, II]) for 
more details. 

• Set-up of the problem. Let T = T n+1 , V be a concavex T-equivariant bundle 
over CP™ of splitting type • • • ; fei, • • • ), and d be a class in 7T2(CP",Z) = Z. Let 
.Mo,o(CP n , cf) be the moduli space of stable maps into CP n , of degree d, genus 0, without 
marked points; and similarly for A7 ,i(CP n , d) and M d = MofiiCP 1 x CP n , (l,d)). Let 
N d = CP( n+1 ) d + n be the linear a-model for CP n . (We will say more about in the next 
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item.) Then one has the following S 1 x T-equivariant diagram: 

V d = 7T*U d U d p*U d V 

J- J- J* •!• 

N d M d M , (CP n ,d) J- A? 0) i(CP n ,d) CP™ 

M^CP 1 xCP",(l,d)) 

where p forgets and eu evaluates at the marked point of a 1-pointed stable map; U d = 
p\ ev* V , the pushforward via p of the pullback of V via ev; it, the contracting morphism, 
is induced by the projection of a stable map in CP 1 x CP™ to the CP™ component and 
contracting the unstable components; and (p, the collapsing morphism, will be explained 
in the next two items. 

Let c top be the top Chern class (i.e. the Euler class) of U d , then the intersection number 
of degree d is defined to be 

K d = I c top (U d ) . 

iM ,o(CP",d) 

One of the goals in the mirror symmetry literatures is to compute K d s and to relate 
them to enumerative problems on CP™, or some varieties therein. An important insight 
from [L-L-Y] is that one can reduce this problem to an easy problem on projective spaces 
{Nd}dLo, cane d t ne linear cr-model in [L-L-Y]. In fact the intersection numbers K d can 
be recovered from cohomology classes Q = {Qdj^Q, called Euler data, defined on those 
projective spaces. In turn the Euler data can be computed essentially by an elementary 
algorithm, and, sometimes, via an explicit formula. A nonlinear recursion involving graph 
sums was use to compute K4 in the case of 0(5) on CP 4 in [Ko]. 

If b is any multiplicative cohomology classes, then we can apply our algorithm to 
compute the integrals 

K d = f_ b(U d ). 

JM ,o(CP n ,d) 

More generally, suppose a, b, . . . are any multiplicative cohomology classes. Then, for any 
given vector bundle V, we have 

a(V) = a (y)+a 1 (V) + ■ ■ ■ + a r {V) , 

where (ii(V) is the degree i component of a(V). We can homogenize a by writing 

a x {V)=Y J * r - l <V), 

i 

where a; is a formal variable, and the class a x remains multiplicative. Likewise, we have 
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etc. Multiplying them, we get a x (V)b y (V) ■ ■ -, which is a new multiplicative class. Note 
that this is a polynomial in the variables x, y, . . ., with coefficients of the form ai(V)bj (V) ■ ■ 
Moreover, each such product correspond to a unique monomial in x, y, .... Now our algo- 
rithm computes K d £ C[x, y, . . .] for the multiplicative class a x (U d )by(U d ) • • • and, hence, 
computes all coefficients of the form 



ai(U d )bj(U d ),... . 

M 0fi (CP n ,d) 

Note that this number is zero unless the integrand has the right total degree. 

Notation. (Cf. Example 1.4.) We shall adapt the following notations for the rest of 
the article: G = S 1 x T, A = (A , A n ), a = ci(C(l)) € H^iBS 1 ), TZ = Q(A)[a], 
n- 1 = Q(A, a), 1ZH* G ( ■) = #£(•) ® Q [A,a] K, TZ^H*^ ■ ) = H* G ( ■ ) ® Q[A>a] TZ' 1 . 

• The linear a-model {N d }f =0 for CP™. Let 

N d = P(F o (CP 1 ,0 CP i(ci)) n+1 ) C P (n+1)d+ri 

be the space of (n + l)-tuple of homogeneous polynomials of degree d on CP 1 up to an 
overall constant multiple in C. An element in N d can be written as 



[^ZQrWQwf T : ••• : ^ZurW^wf 1 

r r 



where [wq : wi] is the homogeneous coordinates for CP 1 and z; ir £ C. The sequence 
{N d }dLo is called the linear a-model for CP™. 

Let G = S 1 x T n+1 ; then G acts on CP 1 x CP n by 

(t, t , t n ) ■ ([w : wi] , [x : ■ ■ ■ : x n \) = {[tw : 101] , [t x : • • • : t n x n }) . 

This induces a G-action on N d with fixed points 

Pi, r = [0 : ••• : : wfrf : : ••• : 0], 

where the non-zero term appears at the i th position, i = 0, • • • , n, and r = 0, ■ ■ ■ , d. Note 
that the weight for the G-action at T Pi r N d is Aj + ra. 

There are two G-equivariant maps between the N d s, defined as follows: 

/ : i\T d _i N d , [fo f n ] ^ [wifo ■ • ■ ■ ■ wifn] ; and 

_ : N d -> N d , [/o^o.fi) : • • • : / n (^o,^i)] ^ [fo(wi,m) : • ■ ■ : fn{wi,w )]. 
From /, one obtaines a chain of inclusions 

N = CP" -U Ni -U ■ ■ ■ -L> N d , 

whose composition gives a canonical incluion I d : Nq = CP™ — > N d . Let k be the equivari- 
ant hyperplane class in H G (N d ). Then the induced map of - on TZ^ 1 H G (N d ) is generated 
by k = k — da, a = —a, and A« = A«. 
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• The G-equivariant morphism 99 : M d — > AT d . First note that and N d are two 

different compactifications of the space A / (o,o( ( CP n , d) of degree d maps from CP 1 to CP n . 
Precisely, an element in _Mo,o(CP n , d) can be written as 

[fo- ■■■ ■ fn] = [^2 z 0rW r w^ r : •■■ : ^ z nr u>o*u;f ~ r ] 

r r 

with /o, • • • , f n relatively prime. Its embedding in N d is tautological, while its embedding 
in M d is given by 

[fo - ■■■ ■ fn] 1 — ► {[wi :wo],[f : ■■■ : /„]) . 

Via these embeddings, the identity map on .Mo,o(CP n , d) extends to a G-equivariant mor- 
phism (p : M d — > TV^. Explicitly, ip can be described as follows: 

Let (/, C) G M d and 7Ti, 7T2 be the projections of CP 1 x CP n onto its first and 
second factor respectively. Then one can decompose C into Co U C\ U • • • UC S 
with Co fl Cj = Xj for j > such that m o f : Co ^ CP 1 and any other Cj 
is pinched to some tt\ o f(xj) = [aj,bj] G CP 1 under 7Ti o /. Let cZj be the 
degree of 7T2 o / : Cj — > CP ra and [<7o : • • • : cr n ] represent the degree do map 
TT2 o / : C -> CP". Then 

</> : (/,C) 1 — ► [ffffo: ••• :5"7n], where g = ]J (ajW - bjW^ . 

3=1 

• Euler data. 

Definition 2.1 [Euler data]. Given an invertible class Q G iJ£(CP n ) -1 , the localization 
of i^(CP n ), an Q-Euler data is a sequence Q = {QdjdLo °f classes G TZH^Nd) that 
satisfy 

(1) Qo = ^. 

(2) The gluing identity: 

i* Pi (n)i;. r (Q d ) = l*Jo^i; io (Q d _ r ), 

for all d and i = 0, ■ • • , n, r = • ■ ■ d. 

An immediate consequence is the following lemma: 

Fact 2.2 [reciprocity]. (Lemma 2.4 in [L-L-Y, I].) // Q is an Euler data, then, for 
i, j = 0, ■ • • , n, r = 0, • ■ ■ , d, d = 0, 1, 2, • ■ one has 

(1) Qd(K + da) = QdJXij. 

(2) Qd{^i)\ a =(\^\ 3 )/d = Qd{^j)\a={\j-Xi)/d for d> 0. 

(3) n{\i) Q d (\j) = Q r {\j)Qd-r{\) ata = (Xj - A;)/r for r > 0. 
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Recall the various bundles and maps from Item Set-up. 



Fact 2.3 [Euler data]. (Theorem 2.8 in [L-L-Y, I].) Let V = V + ®V~ be a concavex 
bundle over CP ra , Xd be the equivariant Euler class ofVd, Qo = Q V = £T{V + )/e T (V~)> 
Qd = <P\{Xd) f or d > 0. Then Q = {Qd} is an Q v -Euler data. 

We call a concavex bundle V — > CP n critical if the induced bundle Ud — > A / io,o(CP ri , d) 
has rank equal to dim (.Mo,o(<CP n , d)) = (n + l)d + n — 3. 

Fact 2.4. (Theorem 3.2 (ii) in [L-L-Y, I].) Let V be a critical concavex bundle over CP". 
Then in the non- equivariant limit A — > 0, 



/ 

J CP" 



-Ht/a hm A ^o J^(Q d ) 
rnU(#-ma)" +1 



a" 3 (2 - <ft) K d . 



Thus, once is determined, the intersection number Kd is also determined. 
Remark 2.5. Since H n+1 = 0, one can rewrite the above formula as 



CP" 



E 

Lk=0 



(-Ht/a) 



k\ 



lim/ d *(Q d ) 



(_l)(n+l)<* 
(rf!)n+lQ,(n+l)d 



E 

fe=0 



(- 



n+l 



a' 3 (2 - efe) if d . 



Note also that in [C-K], there is another formula, implicitly in [L-L-Y, I], that relates Qd 
and Kd,: 



/ 



• Determination of an Euler data. By the localization formula, Qd is determined by 
its restriction i* ir (Qd) at the fixed points pi >r , for i = 0, • • • , n, r = 0, • • • , d. Explicitly, 

Qd = Cp ^ Qd ^ n Q»^,r)(^-Aj -sa) 
(i,r) U(j,s)^(i,r) ( Ai - Aj + (r - a)a ) 

Since the Euler data condition says that 



it turns out that, to determine Q d , one only needs to know its restrictions ip. Q (Qd) at Pi,o 
for i = 0, ■ ■ ■ , n. 

We can now state theorems from [L-L-Y, I] that enables one to determine i*. (Q d )- 

Fact 2.6 [degree bound and determination of Euler data]. (Theorem 2.10, Theorem 
2.11, and Theorem 3.2 (i) in [L-L-Y, I].) Let V be a concavex bundle over CP™ of splitting 
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type (li, I2, ■ ■ ■ ; ki, k 2 , • ■ ■), x be its Euler characteristic class, and Q = {Qd}deNu{o} be the 
X-Euler data for V , as in Fact 2.3. Then the restrictions I^{Qd) £ Hq(CP h ) has 

deg a I* d (Q d ) < (n+l)d-2. 

Furthermore, Q is completely determined by the value of the restrictions ip i0 (Qd)> i = 
0, ■ • • , n, d = 0, 1, 2, • ■ ■, at a = (Xi — Xj)/d> i 3- These values are given explicitly by 



n n ( i * x i 

a m=0 



A,; — A, 



m ■ 



k a d—l 

n n 

b m=l 



m • 



Xi - Xj 



Remark 2.7 [total degree bound]. For V critical, since the rank of U d —> Mo i o('CP n ,d) is 
(n + l)d + n — 3, the total degree of Qd, as a polynomial of k, a, and A, is bounded by 
(n + l)d + n-3. 

These facts and remarks allow one to compute Qd as a polynomial of k and a with 
coefficients in C(Aq, • • • , A n ). We now turn to this detail. 



3 Computation of Qd inductively. 

Following previous notations, let V be a critical concavex bundle over CP n of splitting 
type (l±, I2, • • • ; fci, k 2 , ■••)■ Then Q\ can be computed, using the Atiyah-Bott formula. 
The higher Qd can be computed by the recursive relation from the gluing identities, the 
special values of Qd at the fixed points, and the a degree bound of I d (Qd). 



The computation of Q\. 

For d = 1, deg a Qi(Xi, a) < n — 1 and, for a fixed i, the n-many values Qi(Aj,Aj — Xj), 
j 7^ i are known from Fact 2.6: 

_ la feft-l 

Qi{Xi, Xi - Xj) = n n ( iaXi ~ m {h- Xj)) yi n + m - ^)) > f ° r j ^ * ■ 

a m=0 b m=l 

Thus, using the Lagrange interpolation formula, one obtains 

i* Pi0 (Qi) = Qi(\i,a) = Y, Qi(A,A t -A J ) Uk * Aa ~ Xl+ x \ k) ■ 

3 = 0,-, n llk&jKAk-Xj) 

3 + i 



By the Reciprocity Lemma, i* Pil (Q\) = Q\(Xi,a) = Qi(Aj,— a). In this way, the re- 
striction of Qi at the set of fixed points of the S* 1 x T n+1 -action on N\ are all acquired. 
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Using the localization formula and playing around with the indices, one obtains an exact 
expression 



oi = E 



i=0 



where 



/*(«) 



/i(a) (k - A» - a) + fi(-a) (k - A*) ) JJ(k - A,) J|(/t - A, - a) 



rio-, S )^( 4 ,o)( A * - x j - sa ) 



E 



IL m=o ( - m(A, - Xj) ) n 6 II%=i ( " Mi + m(A t - Xj) ) 
a (a - A, + Aj) IL^(Ai - A*) Ilfc^-^j - x k) 



The computation of for d > 1. 

Let AT = (n + l)eZ + n — 3. Then one may write as a polynomial in k, a with coefficients 
in C[A]: 

N N—fj, 

Qd = E E "V witn "V ^ C[A] . 

M =0 i/=0 

Since we have an explicit formula for Qi, we may assume that Qo, Q\ • • • , Qd-i are all 
determined. 

• Systems from the gluing identity : The gluing identity that an Euler data must satisfy is 
completely encoded in the following two systems of linear equations in w^ v : 

(1) Qd(^i + fa, a) for i = 0, • • • , n and r = 1, • • • , d — 1 : 
The gluing identity says 

Q d (X l +ra, a) = E E "V atl ( x i+ ra T = i P ,,r(Qd) = Qr(K, a) Qd-r{h, a) j Q(A;) 

for i = 0, ■ • • , n and r = 1, ■ • • , d— 1. Denote Q r (Aj, a) Q rf _ r (Aj,a) / J2(Aj) by 
6i(i,r); then &i(i,r) is known by induction, Furthermore, 

deg a (bi{i,r)) < (n + l)r - 2 + (n + l)(d - r) - 2 = (n + l)d - 4 < TV. 

Thus, one may write b\(i,r) = J2f=o b\(i,r, s)a s , where b±(i,r, s) = for 
s > (n + l)d — 3. After expanding the powers and exchanging and relabelling 
the indices to the above equation, one obtains the following linear system in 
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s N-n . . 

E E ^- A r " Ar ~ s = 6i( '' r ' s) 

for i = 0, • • • , n, r = 1, • ■ ■ , d — 1, and s = 0, • • • , N. 



(2) Q d (\i + da, a) for i = 0, • • • , n: 
The gluing identity says 

N N-fj, N N-fj, 

Q d (\i+da,a) = E E "V ^{K+daf = Q d {\i,a) = Q d (X t ,-a) = E E w ^{~ a YK 

fj,=0 v=0 (1=0 v=Q 

for i = 1 = 0, • • • , n. This gives rise to the second linear system in w^ v : 

s-1 N—n . . jV-s , . 

E E /J^Ar- + E^-(i + (-i) s+1 )Ar = o 

for i = 0, • • • , n and s = 0, • ■ ■ , N. 

' System from the special values : 

(3) Q d (\i,^) fori = 0, ■••,n: 

Fact 2.6 gives rise to the third linear system in w;^ : 

e e ^-(^) ^ = nn(w^) 

for i,j = 0, • • • , n, i / j. 

i System from the a-degree bound: 

(4) deg a I* d (Qd) < (n + l)d-2 : 

Since fl£(CP") = C[A, a][«] / n™=o(« " Wd) is obtained by Q d modulo 

the relation Iir=o( K ~~ ^») = 0- This is achieved by iterations of the set of 
replacements 

{K n+l+i ^ K i { _ K n+l + U n =o(K _ Xi)) | i = ,...,JV-n-l} 
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until the K-degree of the resulting Qd is less than n + 1. This can be easily 
done by computer. In this way, one obtains I^(Qd)- Let 



i*AQ d ) = T,< 



2 j Ot t\r • 



Then w'^ is a linear combination of w^ v with coefficients in C[A]. The a-degree 
bound then gives us the fourth system of linear equations in w^ v : 

| w \. = 0, for i > (n + l)d-l. 



Remark 3. 1 . Note that the whole content of gluing identities in the definition of Euler data 
is already absorbed in the first and the second linear systems (1) and (2) above. Since all 
the identities that appear in the Reciprocity Lemma are obtained by substituting into the 
gluing identities some special a values, they will be automatically satisfied once System 
(1) and System (2) above are satisfied. Thus, they do not provide us with extra equations. 
Furthermore, Fact 2.6 implies that the above system has a unique solution. 



4 Modifications for non-critical bundles over CP n . 

So far, our discussion has been focusing on critical bundles. For non-critical bundles V 
over CP™, rank(Ud) and dim Alo,o(CP n , d) are different for some d, by definition. If one 
takes {Qd}d associated to the top Chern class of V, then one will simply get for Kd- Let 
r be the rank of V. Then, to obtain more interesting invariants, one may consider taking 
{Qd}d to be the Euler data associated to the Chern polynomial 

c( x ) = x r + c 1 {V)x r - 1 + +c top (V) 

of V. In this case, some details in Sec. 2 and Sec. 3 will have to be modified accordingly; 
however, the conceptual flow remains the same. 

The various items that are involved in the actual computation of Qd and Kd and their 
modifications following [L-L-Y] are listed below: 

• degree bound of Qd ■ Since the top Chern class gives the highest degree terms, 
both the total degree bound and the a-degree bound remain valid. 

• the special values iL.AQd) a,-a, : These special values are now given by 

f»,u a= _J_ 

^(Qd)^ = nil {x + lah-m^^)\{ b Y{ (x-k b \ t + m^^). 

d a m=0 b m=l 

This leads to a corresponding change in the computation of Q\ and right-hand-side 
of the third linear systems in Sec. 3. The first, second, and the fourth linear systems 
remains valid. 
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• from Q d to K d : Let s = rank(U d ) - dim M(CP n , d) . Theorem 7.2 in [L-L-Y, II] 

gives 

- — f c~ m/a lim A^*(Qrf) = J_ (2 _ dt)K . 

si dx s x=0 7cp™ nl=i (# - « 3 ^ s 



Note that, for the validity of this formula alone, it is not required that s d be indepen- 
dent of d [Li]. Thus, once Q d is obtained by solving the system of linear equation, 
K d follows. 

Based on Sec. 3 and the discussion here, a code EulerData_MP .m is written. The detail 
is in Sec. 6. 



5 Examples. 



Using the Maple code "EulerData_MP .m" in Sec. 6 and taking 0, to be the Chern polyno- 
mial, we compute the first few K d for some non-critical bundles, as listed in Table 5-1 
(cf. Cases 6, 14-18 in SEC. 3 of the code). 



bundle 


d = 1 


d = 2 


d = 3 


T,CP 2 


10 x 2 






C p 4 (6) 


50400 x 


(752729895/4) x 2 


(433244745198080/243) x 3 


O cp4 (7) 


451570 x 2 


(403985396325/32) x 4 


(15755269694706695755/17496) x 6 


O CP 4(8) 


2773820 x 3 


(3178734062035/8) x 6 


(46028387589557254161275/314928) x 9 


cp4 (9) 


13198850 x 4 


(243281907041715/32) x 8 


(197802281929974511821535/17496) x 12 


O cp4 (10) 


52040450 x 5 


(25908993204089625/256) x 10 


(71418501571607082433686025/139968) x 15 



Table 5-1. K d for some non-critical bundles. 



Note that, in Table 5-1, if one formally converts K d to n d using the usual multiple 
cover formula K d = J2k\d then the n d thus obtained will no longer be integers for 

d>2. Compared with the other twelve examples computed/tested using the code, this in- 
dicates that for those non-critical bundles over CP™, whose rank (U d ) — dim (A / io,o(CP n , d)) 
depends on d in a non-trivial way, the above multiple cover formula will have to be mod- 
ified. Exactly how is an issue for further investigation. 

6 A package in Maple V for the computation of Qd and K^. 

General remarks. 

A few remarks are given below concerning the code, its current scope, and its usage. 
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• Two bad things first. The first one is that the code, as currently written, is limited 
only to the case that Q is an Euler class or a Chern polynomial. For other multi- 
plicative characteristic classes, one will not be able to use the code to make sensible 
computation without entering the core part, SEC. 1.3 (Q1K1) and SEC. 2 (Euler- 
Data), to do some modifications. This shortage will hopefully be gradually removed 
along with the development of the theory, using the relation of a given multiplicative 
class with Chern roots. The second one is that the actual computation of Euler data 
and could be a very demanding task both for Maple V and the machine (cf. some 
words in SEC. 3 of the code about Test 6: T^CP 2 even for d = 2). When it exceeds 
the capacity the Maple V, one will likely get an error message. Experience tells us 
that occasionally these messages may be misleading. 

• Now let us mention something more positive. As indicated in SEC. 3 of the code, the 
code has been tested correct for all known cases within the capacity of Maple V and 
the machine used. It is also made user- friendly: to run for a case of study, one only 
has to follow the examples in SEC. 3 and modify the various arguments/parameters 
to be fed into the function 'EulerData', as instructed there. By no means does one 
need to do anything else; nor do we assume any knowledge of Maple at all. 

• The specialization used, the time consumed, and the RAM memory used for the 
examples tested, particularly those that take long hours, are recorded in SEC. 3 of 
the code for reference. Only the last run of each case is recorded in these notes. 



Instruction of running the code under Window 98. 

For non-maple-user, let us give here some instruction of running the code under Window 
98: (assuming there is already Maple V in the folder "Programm Files") 

• Read first the instructions both at the beginning of the code and at the start of SEC. 
3 of the code . 

• Save the code (say by the file name 'EulerData_MP .m') in the subfolder 'Bin. wnt' in 
the folder 'Maple V Release 5' that is automatically set up in the folder 'Programm 
Files' when installing Maple V. Double click the maple icon to open a Maple V work- 
sheet. Inside the worksheet, type in the following command line after the prompt 

restart; read "EulerData_MP .m" ; 



then hit Enter on the keyboard. The output will be displayed directly on the 



worksheet. 

For all other operating systems, please consult the system manager. 
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The code 'EulerData_MP .m\ 



#========================( Beginning of the code. )======================== 

# ( PLEASE READ ME BEFORE USE ) 

# NOTE 1. Though the code has been thoroughly tested, nevertheless if you find any bug that has escaped 

# our scrutiny, please send an e-mail to 'chienliuOmath.harvard.edu' for remedy of this. 
# 

# NOTE 2. CAUTION TO USER: 

# The machine used for the tests is 
# 

# Dell PC with Pentium II 400 MH, Cache 512 KB, RAM 384 MB, Hard Drive 7 GB. 
# 

# Make sure the code is run in a machine at least of the same level as above to avoid possible crashes. 
# 

# NOTE 3. Anything free comes with no guarantee. Please take your own caution/ judgement or consult experts 

# BEFORE running the code. The author of the code shall not be responsible for any possible damage 

# and/or loss caused by the code. The author does not recommend the use of the code for applications 

# in which errors or omissions could threaten life, injury or significant loss. 
#========================================================================== 



# 

# Language used: Maple V, created by Waterloo Maple Inc. at Waterloo, Ontario, Canada. 
# 

# 

# SEC. Introduction and outline. 

# 

# 

# EulerData_MP.m — a package in Maple V. 

# Given an Euler class or the Chern polynomial of a concavex bundle $V$ over $CP~n$, the 

# package computes the Euler data $Q=\{Q_d\}$ and the intersection numbers $K_d$ for $V$ 
# 

# Date of completion: November 21, 1999. 

# Test: All the subroutines involved are tested correct separately. See the individual remarks for 

# details. Some records of the run are recorded in SEC. 3 for a reference of future improvements. 

# Date of last revision: January 7, 2000. 

# 

# List of functions/subroutines: 

# Main routine (Sec. 2): EulerData. 

# Subroutines (Sec. 1): LSolveSub (Sec. 1.1.1), LSolveSub_A (Sec. 1.1.2 ), 

# PolyCongruent (Sec . 1 . 2) , Q1K1 (Sec. 1.3). 

# 

# SEC. 3 contains INSTRUCTIONS of using the code and various cases of computation and serves as 

# a self-explanatory TUTORIAL for users. 

# NOTE. For the computataion of $K_d$ alone, one may add a specialization BEFORE employing the function 

# 'EulerData' . Different cases or different degrees for the same case may require different 

# specializations. To be compatible with the code, the specialization must be of the form in the 

# variables 'lambda'and 'u' . The specialization is given in the form: 
# 

# lambda :i -> some_made_up_f unction(i , u) 
# 

# such that lambda(i)=0 when u=0. (See the examples in SEC. 3.) 

# 



# 

# SEC. 1 Subroutines. 
# 
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# 

# SEC. 1.1.1 Definition of the function 'LSolveSub' . 

# 

# ' LSolveSub (equation_set , variable_list , check)' 

# solves a system 1 equation_set ' of linear equations in variables in 'variable_list' by 

# eliminations and back substitutions. 
# 

# NOTE: If 1 check' =0 (resp. 1), then the consistency of the system will not (resp. will) be checked. 
# 

# Date of completion: November 19, 1999. 

# Test: Tested correct for six sample examples. 

# Date of last revision: November 23, 1999. 

# 

LSolveSub := proc( equation_set , variable_list , check ) 
local eqn_list_diag, eqn_list_tri , eqn_set, i, j, list_part, 
m, modification, s, sol, t, test, v, vl: 

eqn_set := equation_set : 

eqn_set := map( el->lhs(el)-rhs(el)=0, eqn_set ): 
m := nops(eqn_set) : 
v := variable_list : 
vl : = nops (v) : 
eqn_list_tri := [] : 

# print("Number of equations = ", m) : 

# print("Number of variables = ", vl) : 

# print ("LSolveSub: 1000000"): 

for i from 1 to vl do 

# printC'i = ", i) : 

for j from 1 to m do 

# print("j = ", j): 
t := j: 

test := coeff( lhs(eqn_set [j] ) , v[i] ): 

# print ("v[i] = ", v[i]): ################# 

# print ("LSolveSub: 1220000"): ################# 
if testoO then 

# printC'test = ", test): ################# 

# print ("eqn_set[j] = ", eqn_set [j] ) : ################# 

# print ("v[i] = ", v[i]) : 

s := solve ( { eqn_set[j] }, { v[i] } ): 

# print ("s = " , s) : 
sol := s[l] : 

f i: 

# print ("LSolveSub: 1240000"): ################# 
if testoO then break fi: 

od: 

# print ("LSolveSub: 1400000"): 

eqn_list_tri := [ op(eqn_list_tri) , sol ]: 
eqn_set := eqn_set minus { eqn_set [t] }: 

if nops(eqn_set)=0 then break fi: 

eqn_set := map( el -> subs( sol, el ), eqn_set ): 
eqn_set := map( el -> expand( lhs(el) )=0, eqn_set) : 
eqn_set := eqn_set minus {0=0}: 
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if eqn_set <> {} then 

eqn_set := map( el -> collect( lhs(el), v )=0, eqn_set ) 

fi: 

m := nops( eqn_set ): 
od: 



# 

# In terms of matrices, so far an upper triangular matrix is formed. 

# The next step is the back substitutions to solve v. 

# 

# print ("LSolveSub: 6000000"): 

m := nops( eqn_list_tri ): 
eqn_list_diag := eqn_list_tri : 

unassign( ' eqn_list_tri ' ): 

for i from 2 to m do 

list_part := [ op(- (i-1) . . -1 , eqn_list_diag)] : 
modification := rhs( op( -i, eqn_list_diag ) ): 
modification := subs( list_part, modification ): 
modification := normal( modification ): 
modification := simplify( modification ): 

eqn_list_diag := subsop( [-i , 2] modification, eqn_list_diag ): 
od: 

# 

# consistency check. 

# 

# print ("LSolveSub: 9000000"); 

if check=l then 

if eqn_set<>{} then 

eqn_set = map( el -> subs( eqn_list_diag, el), eqn_set ): 
eqn_set = { map( el -> lhs(el)-rhs(el) , eqn_set ) }: 
if eqn_set <> {0} then 
print ( " inconsistent " ) 
else print ("consistent") 

fi: 

else print("consistent") 

f i: 

f i: 

unassign( ' eqn_list_tri ' , 'eqn_set', 'i', 'j', 'list_part', 

'm', 'modification', 's', 'sol', 't', 'test', 'v', 'vl'): 
RETURN (eqn_list_diag) : 
end: 



# 

# SEC. 1.1.2 Definition of the function 'LSolveSub_A' . 

# 

# 'LSolveSub_A(equation_set, variable_list , check)' 

# solves a system ' equation_set ' of linear equations in variables in 'variable_list' 

# by eliminations and back substitutions. 
# 

# NOTE: If ' check '=0 (resp. 1), then the consistency of the system will not (resp. will) be checked. 
# 
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# Date of completion: November 19, 1999. 

# Test: Tested correct for six sample examples. 

# Date of last revision: January 6, 1999. 

# Modification from 'LSolveSub' : 

# Replacing the built-in function 'solve' by direct algebra to increase the capacity of the code. 
# 

LSolveSub_A := proc( equation_set , variable_list , check ) 
local eqn_list_diag, eqn_list_tri , eqn_set, i, j, list_part, 
m, modification, s, s_l, sol, t, test, v, vl: 

eqn_set := equation_set : 

eqn_set := map( el->lhs(el)-rhs(el)=0, eqn_set ): 
m := nops(eqn_set) : 
v := variable_list : 
vl := nops(v) : 
eqn_list_tri := [] : 

# print ( "Number of equations = ", m ): 

# print ( "Number of variables = ", vl ): 

# print( "LSolveSub_A: 1000000" ): 

for i from 1 to vl do 

# print( "i = ", i ): 

for j from 1 to m do 

# print ( "j = ", j ): ##################### 

# ( "LSolveSub_A: 1220000" ): ########## 
t := j: 

test := coeff( lhs(eqn_set [j] ) , v[i] ): 

# print ( "v[i] = " , v[i] ) : 
if testoO then 

# print ( "test = ", test ): ##################### 

# 

# Remark: Replacement of the following line in 'LSolveSub'. 

# 

# s := solve ( { eqn_set[j] }, { v[i] } ): 

# 

# print ( "v[i] = ", v[i] ): 

s_l := -(lhs(eqn_set [j] ) -test*v [i] ) /test : 
s_l := normal (s_l): 
s := {v[i]=s_l}: 

# print ( "s = ", s ) : 

# 

# End of replacement . 

# 

sol := s[l] : 

f i: 

if testoO then break fi: 

od: 

# print( "LSolveSub_A: 1400000" ): 

eqn_list_tri := [ op(eqn_list_tri) , sol ]: 
eqn_set := eqn_set minus { eqn_set [t] }: 

if nops(eqn_set)=0 then break fi: 

eqn_set := map( el -> subs( sol, el ), eqn_set ): 
eqn_set := map( el -> expand( lhs(el) )=0, eqn_set) : 
eqn_set := eqn_set minus {0=0}: 



18 



if eqn_set <> {} then 

eqn_set := map( el -> collect( lhs(el), v )=0, eqn_set ) 

fi: 

m := nops( eqn_set ): 
od: 

# 

# In terms of matrices, so far an upper triangular matrix is formed. 

# The next step is the back substitutions to solve v. 

# 

# print ( "LSolveSub_A: 6000000" ): 

m := nops( eqn_list_tri ): 
eqn_list_diag := eqn_list_tri : 

for i from 2 to m do 

list_part := [ op(- (i-1) . . -1 , eqn_list_diag)] : 
modification := rhs( op( -i, eqn_list_diag ) ): 
modification := subs( list_part, modification ): 
modification := normaK modification ): 
modification := simplify( modification ): 

eqn_list_diag := subsop( [-i , 2] modification, eqn_list_diag ): 
od: 

# 

# consistency check. 

# 

# print ("LSolveSub_A: 9000000"): 

if check=l then 

if eqn_set<>{} then 

eqn_set = map( el -> subs( eqn_list_diag, el), eqn_set ): 
eqn_set = { map( el -> lhs(el)-rhs(el) , eqn_set ) }: 
if eqn_set <> {0} then 
print ( " inconsistent " ) 
else print ("consistent") 

fi: 

else print ("consistent") 

fi: 

f i: 

unassign( ' eqn_list_tri ' , 'eqn_set', 'i', 'j', 'list_part', 

'm', 'modification', 's', 's_l', 'sol', 't', 'test', 'v', 'vl'): 

RETURN (eqn_list_diag) : 
end: 



# 

# SEC. 1.2 Definition of the function 'PolyCongruent ' . 

# 

# 'PolyCongruent (polynomial, replacememt , variable)' 

# gives the polynomial of minimal degree in 'variable' that is congruent to 'polynomial' 

# modulo the relation 'replacement'. 
# 

# Note: The argument 'replacement' must be of the form: 
# 

# variable '"n = polynomial in 'variable' of degree <= n 
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# 

# 

# Date of completion: November 20, 1999. 

# Test: Tested correct for sample examples. 

# Date of last revision: November 20, 1999. 

# 

PolyCongruent : = proc ( polynomial , replacement , variable ) 
local dO, dl, dnow, i, relation, relation_set , poly, v: 

poly := polynomial: 
relation := replacement: 
v := variable: 

dO := degree( Ins (relation) , v ): 
dl := degree ( poly, v ): 

relation_set := 

{ seq( v~(i+dO)=collect( v~i*rhs (relation) , v ), i=0..dl-d0 ) }: 
dnow : = dl : 

for i from 1 to dl-dO+1 do 

if dnow < dO then break fi: 

poly := subs( relation_set , poly ): 

poly := expand ( poly ): 

poly := simplify( poly ): 

poly := collect( poly, v ): 

dnow := degree ( poly, v ): 
od: 

unassign( 'dO' , 'dl', 'dnow', 'i', 'relation', 'relation_set' , 'v' ): 

RETURN (poly) : 

end: 



# 

# SEC. 1.3 Definition of the function 'Q1K1'. 

# 

# 'QlKl(n, splitting_type, s_diff_0, opt_l) ' 

# gives the first term in the $0mega$-Euler data and the first Kontsevich number for a concavex 

# bundle $V$ over ${\Bbb C}P~n$, where $0mega$ is the Chern polynomial. 
# 

# Note.l: The argument ' splitting_type ' must be of the form "[[_, _, ...],[_, _, ...]]", where the 

# first [_, _, ...] comes from the convex summand and the second [_, _, ...] comes from the 

# concave summand of the concavex bundle $V$. 
# 

# Note. 2: The argument 'opt_l' takes only values and 1. 

# If 'opt_l' = 0, then only $K_1$ will be computed. 

# If 'opt_l' = 1, then exact $Q_1$ will also be computed. 

# 

# Date of completion: November 19, 1999. 

# Test: The value Kl are correct for nine known cases: Cases 1-5, 7-10 in SEC. 3. 

# Date of last revision: November 19, 1999. 

# 

Q1K1 := proc(n, splitting_type , s_diff_0, opt_l) 
local i, instl, k, kl, kla, 11, 12, poly, polyl, 

qlpolyO, qlpolyOl, qlpoly02, qlpolylO, qlpolyll, qlpolyllO, qlpolylll, 

qlsummandl, qlsummand2, qlvalue, qlvaluel, qlvalue2, 

ql, qlO, qlOa, qll, qlla, si, s2, v, y: 

11 := splitting_type [1] : 

12 := splitting_type [2] : 
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si := nops(ll) : 
s2 := nops(12) : 

printC'll = ", 11) : 

print ("12 = ", 12) : 

print("sl = " , si) : 

print ("s2 = " , s2) : 

printC'QlKl, 1000000"): 



qlvaluel := proc(v, y) 
local a, m: 
if sloO then 

mul( mul(x+ll[a]*v-m*y, m=0 . . 11 [a] ) , a=l..sl ) 
else 1 

fi: 
end: 

qlvalue2 := proc(v, y) 
local b, m: 
if s2<>0 then 

mul( mul(x+12[b]*v+m*y, m=l . . -12 [b] -1) , b=l..s2 ) 
else 1 

fi: 
end: 

qlvalue := (v, y) -> qlvaluel (v, y)*qlvalue2(v, y) : 

qlsummandl := proc(i, j, alpha) 
local factorl, factor2, k: 
factorl := 1: 
for k from to n do 

if koi and koj then 

factorl := factorl*( alpha-lambda(i)+lambda(k) ) 
else factorl := factorl 

fi: 

od: 

factor2 := 1: 

for k from to n do 

if koi and koj then 

factor2 := factor2*( lambda (k) -lambda(j ) ) 

else 1 

fi: 

od: 

qlvalue( lambda(i) , lambda(i)-lambda(j ) ) *f actorl/f actor2 : 
end: 

qlsummand2 := (i, j, alpha) -> 

normal ( qlsummandl (i , j, alpha) ): 



# 

# THEORY: 'qlpolyOKi, alpha)' = $i~{\ast}_{p_{i,0}}(Q_l)$ 
# 

qlpolyOl := proc(i, alpha) 
local j, t: 
t := 0: 

for j from to n do 
if j<>i then 

t := t + qlsummand2(i , j, alpha) 
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else t := t 

f i: 
od: 
end: 

qlpoly02 := (i, alpha) -> normaK qlpolyOKi, alpha) ): 
qlpolyO := (i, alpha) -> simplify( qlpoly02(i, alpha) ): 

qlpolylO := (i, alpha) -> normaK qlpolyO(i, alpha) ): 
qlpolyllO := (i, alpha) -> simplify( qlpolylO (i, alpha) ): 
qlpolyll := (i, alpha) -> qlpoly 10 (i, -alpha) : 
qlpolylll := (i, alpha) -> simplify( qlpolyll (i , alpha) ): 

# 

# definition of qlO 

# 

qlOa := proc(n) 

local i, ql0_factor2, ql0_factor3, total: 

ql0_factor2 := proc(i, n) 
local j , k, t : 
t := 1: 

for j from to n do 
for k from to 1 do 
if joi or k<>0 then 

t := t*(kappa-lambda(j)-k*alpha) 
else t:= t 

f i: 
od: 
od: 
end: 

ql0_factor3 := proc(i, n) 
local j , k, t : 
t := 1: 

for j from to n do 
for k from to 1 do 
if joi or k<>0 then 

t := t*( lambda(i)-lambda(j)-k*alpha ) 
else t := t 

f i: 
od: 
od: 
end: 

total := add( qlpolyll0(i,alpha)*ql0_f actor2(i,n)/ql0_f actor3(i,n) , 
i=0. .n ) : 

end: 

qlO := normaK qlOa(n) ): 
qlO := simplify( qlO ): 

# 

# defintion of qll 

# 

printC'QIKI, 4000000"): 

qlla := proc(n) 

local i, qll_factor2, qll_factor3, total: 
qll_factor2 := proc(i, n) 
local j , k, t : 
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t := 1: 

for j from to n do 
for k from to 1 do 
if joi or kol then 

t := t*( kappa-lambda (j ) -k*alpha ) 
else t := t 

f i: 

od: 
od: 
end: 

qll_factor3 := proc(i, n) 
local j , k, t : 
t := 1: 

for j from to n do 
for k from to 1 do 
if joi or kol then 

t := t*( lambda(l)-lambda(j)+(l-k)*alpha ) 
else t := t 

f i: 
od: 
od: 
end: 

total := add( qlpolylll (i , alpha) *qll_f actor2(i, n)/qll_f actor3(i,n) , 
i=0. .n ) : 

end: 

qll := normaK qlla(n) ): 
qll := simplify( qll ): 

printC'QlKl, ") : 

# 

# The first term of Euler data : Ql. 

# 

ql := qlO+qll: 
ql := normal (ql) : 
ql := simplify (ql) : 

printC'QlKl, 8000000"): 

# 

# Computation of Kl from Ql, using the formula in Lian, Liu, and Yau's "Mirror Principle I, II" 

# and [Cox-Katz] "Mirror symmetry and algebraic geometry" . 

# 

poly := ql: 

if opt_l=0 then 

poly := subs( u=0, ql ) 

else 

for i from to n do 

poly := subs( lambda(i)=0, poly ) 
od: 

fi: 

poly := subs( kappa=h, poly ): 

# 

# Extra manipulation for noncritical bundles with 0mega=the Chern polynomial. 
# 



23 



if s_diff_0(l) <> then 

poly := diff( poly, x$s_dif f _0(1) ): #( Cf . Mirror Principle II, p. 36. ) 
poly := subs( x=0, poly ) 

f i: 

# 

polyl := add( (-h*t/alpha) ~k/(k ! ) , k=0..n )* (-1) " (n+1) /(alpha" (n+1) ) * 

(add( (h/alpha)"k, k=0 . .n) ) ~ (n+1) *poly : 
polyl := simplify (polyl) : 

kl := alpha~3/(2-t)*coeff (polyl, h~n) : 
kl : = normal (kl) : 
instl := kl: 

# 

# Extra manipulation with Omega=the Chern polynomial 

# 

if s_diff_0(l) <> then 

instl := kl/(s_diff_0(l) !) : 

kl := x-(s_diff_0(l))*kl/(s_diff_0(l) !) 

fi: 

print ( "Routine : Q1K1") : 
print ("kl = ", kl) : 
print ("instl = ", instl): 

# 

# kla := alpha~2*coeff (h*polyl, h~n) : 

# kla := normal (kla): 

# 

unassign( ' i ' , 'k', '11', '12', 'poly', 'polyl', 'qlpolyO', 'qlpolyOl', 
'qlpoly02', 'qlpolylO' , 'qlpolyll', 'qlpolyllO' , 'qlpolylll', 
'qlsummandl ' , 'qlsummand2' , 'qlvalue', 'qlvaluel', 'qlvalue2' , 
'qlO', 'qlOa', 'qll', 'qlla' , 'si', 's2', 'v' , 'y' ): 

RETURN ( [ql, kl , instl] ): 
# RETURN ([ql, kl , kla]): 
end: 



# 

# SEC. 2 The MAIN ROUTINE. 

# 

# 

# Definition of the function 'EulerData' . 

# 

# 'EulerData(n, splitting_type , Omega_0, s_diff_0, dmax, opt_l, opt_2, check)' 

# computes $Q_1 Q_d$ with $Q_0=\0mega_0$ and $K_1 K_d$ for a convex bundle $V$ over $CP"n$ 

# of splitting type ' splitting_type ' . 
# 

# Note.l: The argument ' splitting_type ' must be in the form of lists "[[_, _, ...],[_, _, ...]]", where 

# the first [_, _, ...] comes from the convex summand and the second [_, _, ...] comes from the 

# concave summand of the concavex bundle $V$. 
# 

# Note. 2: The argument 's_diff_0' is the differnce: rank(U_d) -dim(\bar{\cal M}_{0,0}(CP~n, d)); 

# (cf. [L-L-Yl, 2] for notations) . 
# 
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# Note. 3: The argument 'opt_l' takes only values and 1. 

# If 'opt_l' = 0, then only $K_d$, $d=l , .., dmax$, will be computed. 

# If 'opt_l' = 1, then exact $Q_d$ up to $Q_{dmax}$ will also be computed. 
# 

# Note. 4: The argument 'opt_2' takes only values and 1. 

# If 'opt_2' = 0, then 'LSolveSub' is enmployed; 

# if 'opt_2' = 1, then 'LSolveSub_A' is employed. 
# 

# Note. 5: The argument 'check' takes only values and 1. 

# If 'check' =0 (resp. 1), then consistency check of the linear system will not (resp. will) 

# be performed. 

# 

# Date of completion: November 21, 1999. 

# Test: See SEC. for details. 

# Date of last revision: December 8, 1999. 

# 

EulerData := proc(n, splitting_type , 0mega_0, s_diff_0, dmax, opt_l, check) 
local d, eqn_set, inhomogenll, inhomogenl2, inhomogenl3, inhomogenl, 

inhomogen3, inhomogen31, inhomogen32, instanton_list , 

k, kl, kla, ka, kd, kda, kdn, 11, 12, nl, 

poly, polyl, poly_set, poly_setl, poly_setll, poly_set21, 
poly_set2, poly_set3, poly_set4, q, qO, ql , ql_dif f erence , 
qd, qdll, qdlll, qdll2, 

qdll3, qdll4, qdl2, qdl3, qd21, qd22, qd23, qd31 , qd311, qd312, qd313, 
qd32, qdr, qdrl, qd_value, qd_valuel, qd_value2, qkl , qkl_linear, 
rel, si, s2, s_diff, substitution, variable_list : 

qkl := QlKl(n, splitting_type , s_diff_0, opt_l) : 

print("d=l") : 

print ( "kl = " , qkl [2] ) : 

# printC'kl =", qkl [2] , qkl [3] ) : 

print ("x = " , x) : 

print ("Routine: EulerData; Omega = ", Omega): 
if dmax=l then RETURN (qkl) fi: 

# 

# Computations of $Q_d$ and $K_d$ for $d>=2$. 

# 



11 
12 
si 
s2 



= splitting_type [1] : 
= splitting_type [2] : 
= nops (11) : 
= nops (12) : 



qO_inverse := l/0mega_0: 

print ("Routine: EulerData; qO_inverse = ", qO_inverse) : 



ql := qkl[l] : 
q := [ ql ]: 
kl := qkl [2] : 
k := [ kl ] : 

# kla : = qkl [3] : 

# ka := [ kla ] : 
instanton_list := [ qkl [3] ]: 

print("000000000000000") : 
print ("k = ", k) : 

# print ("ka = ", ka) : 
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print ("instanton_list = ", instanton_list) : 



for d from 2 to dmax do 
printC'd = ", d) : 
nl := (n+l)*d+n-3: 

qd := add( add( w(mu,nu) *alpha~mu*kappa~nu, nu=0..nl-mu ), mu=0..nl ): 
variable_list := [ seq( seq( w(mu, nu) , nu=0 . . nl-mu ), mu=0..nl ) ]: 



# 

# Linear system from the 1st set of equations. 

# 

print("MAIN: 1000000; I am now working out the first set of equations."): 



qdlll := (i, r) -> subs( kappa=lambda(i)+r*alpha, 

qdll2 := (i, r) -> expand( qdllld, r) ): 

qdll3 := (i, r) -> simplify( qdll2(i, r) ): 

qdll4 := (i, r) -> collect( qdll3(i,r), alpha ): 



qd ): 



inhomogenll := (i, r) -> 

subs( { kappa=lambda(i) , alpha=-alpha} , q[r] )* 

subs( kappa=lambda(i) , q[d-r] )*subs( h=lambda(i) , q0_inverse ): 



inhomogenl2 
inhomogenl3 
inhomogenl 



(i, r) -> normaK inhomogenll (i , r) ): 

(i, r) -> simplify( inhomogenl2(i , r) ): 

(i, r) -> collect( inhomogenl3(i, r) , alpha ): 



qdll := (i, r) -> qdll4(i, r)-inhomogenl (i , r) : 
qdl2 := (i, r) -> simplify( qdll(i, r) ): 
qdl3 := (i, r) -> collect( qdl2(i, r) , alpha ): 
poly_setll := (i, r) -> coeffs( qdl3(i, r) , alpha ) : 
poly_setl := { seq( seq( poly_setll (i , r) , r=l..d-l ), 
poly_setl 



i=0. .n ) }: 

:= map( el -> collect( el, variable_list ), poly_setl ): 



# 

# Linear system from the 2nd set of equations. 

# 

print("MAIN: 2000000; I am now working the second set of equations."): 



qd21 := i -> subs( kappa=lambda(i)+d*alpha, qd )- 

subs( { kappa=lambda(i) , alpha=-alpha }, qd ) : 
qd22 := i -> simplify( qd21(i) ): 
qd23 := i -> collect ( qd22(i), alpha ): 
poly_set21 := i -> coeffs( qd23(i), alpha ): 
poly_set2 := { seq( poly_set21 (i) , i=0..n ) }: 

poly_set2 := map( el -> collect(el, variable_list ), poly_set2 ): 

# 

# Linear system from the 3rd set of equations: 

# the special values $Q_d(\lambda_i , (\lambda_i-\lambda_j)/d )$ 

# 

print("MAIN: 3000000; I am now working out the third set of equations."): 

qd_valuel := proc(v, y) 
local a, m, t: 
if slOO then 

t := mul( mul( x+11 [a] *v-m*y , m=0..11[a]*d ), a=l..sl ) 
else t := 1 

fi: 

RETURN (t) : 
end: 



qd_value2 := proc(v, y) 
local b, m, t: 
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if s2<>0 then 

t := mul( mul( x+12 [b] *v+m*y , m=l . . -12 [b] *d-l ), b=l . . s2 ) 
else t := 1 

fi: 

RETURN (t) : 
end: 

qd_value := (v, y) -> qd_valuel(v, y)*qd_value2(v, y) : 

inhomogen31 := (i, j) -> qd_value( lambda(i) , (lambda(i) -lambda( j ) ) /d ): 
inhomogen32 := (i, j) -> normaK inhomogen31 (i , j ) ): 
inhomogen3 := (i, j) -> simplify( inhomogen32 ( i , j) ): 
qd311 := (i, j) -> 

subs( { kappa=lambda(i) , alpha=( lambda(i) -lambda( j ) )/d }, qd) : 
qd312 := (i, j) -> expand( qd311(i, j) ): 
qd313 := (i, j) -> simplify( qd312(i, j) ): 
qd31 := (i, j) -> qd313(l, j)-inhomogen3(i, j): 
qd32 := (i, j) -> simplify ( qd31(i, j) ): 

poly_set3 := { seq( seq( qd32(i, j), j=i+l..n), i=0..n ) }: 
poly_set3 := map( el-> collect( el, variable_list ), poly_set3 ): 

# 

# Linear equations from the 4th set of equations: 

# $\alpha$-degree bound: $\deg_{\alpha} I~{\ast}_d(Q_d) \le (n+l)d-2$. 
# 

print("MAIN: 4000000; I am now working out the fourth set of equations."): 

rel := mul( kappa-lambda(i) , i=0..n ): 
rel := expand( rel )-kappa~ (n+1) : 
rel := collect ( rel, kappa ): 
rel := kappa" (n+l)=rel : 

qdr := PolyCongruent ( qd, rel, kappa ): 
qdr := collect ( qdr, alpha ): 

poly_set4 := { seq( coeff( qdr, alpha~i ), i=(n+l) *d-l . .nl ) }: 
poly_set4 := map( el -> collect( el, kappa ), poly_set4 ): 
poly_set4 := map( el -> coeffs( el, kappa ), poly_set4 ): 

# 

# Combination of the four linear systems and solve the system. 

# 

print("MAIN: 6000000; I am now combining the four systems of equations."): 

poly_set := 'union' ( poly_setl, poly_set2, poly_set3, poly_set4) ; 
eqn_set := map( el -> el=0, poly_set ); 

unassign( 'poly_set ' , 'poly_setl', 'poly_set2', 'poly_set3', 'poly_set4' ) : 

## 

## Added for testing the consistency check routine. 
## 

## print( "op(-l, eqn_set) = ", op(-l, eqn_set) ): 

### eqn_set := { op( eqn_set ), lhs( op(-l, eqn_set) )=100 }: # or 
### eqn_list := { op( eqn_set ), 2*lhs( op(-l, eqn_set) )=0 }: 
## 

## End of testing lines 

## 

if opt_2 = then 

print("I am now solving the system of linear equations using 'LSolveSub' 
substitution := LSolveSub( eqn_set, variable_list , check ) 
else 

print ("I am now solving the system of linear equation using 'LSolveSub_A 
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substitution := LSolveSub_A( eqn_set, variable_list , check ) 

fi: 

print("substitution = ", substitution): 
qd := subs( substitution, qd ); 
qd := normal (qd) ; 
qd := simplify (qd) ; 
q := [ op(q) , qd ] ; 

print("The system of linear equations is now solved."): 
print ("q = " , q) : 

# 

# The nonequivariant limit of $I_d"{\ast} (Q_d) $ . 

# 

print("I am now turning to compute Q_d, K_d, and n_d."): 

qdrl := PolyCongruent ( qd, rel, kappa ): 

if opt_l=0 then poly := subs( u=0, qdrl ) 
else poly := qdrl: 

for i from to n do 

poly := subs( lambda(i)=0, poly ) 

od: 

fi: 

poly := subs( kappa=h, poly ): 

# 

# Extra manipulation for noncritical bundles with Omega=the Chern polynomial. 

# (Cf. Mirror principle, I: p. 37.) 

# 

if s_diff_0(d) <> then 

poly := diff( poly, x$s_dif f _0(d) ): 
poly := subs( x=0, poly ) 

fi: 

# 

# Computation of $K_d$ via $Q_d$, using the formula in 

# [Lian-Liu-Yau] "Mirror principle I", "Mirror principle II". 

# 

polyl := add( (-h*t/alpha) ~r/ (r ! ) , r=0..n )* 

mul( ( add( (h/ (m*alpha) ) ~r , r=0..n ) )~(n+l), m=l..d )* 

poly/( (d!)-(n+l)*(-alpha)-(d*(n+l)) ): 
polyl := normal(polyl) : 
polyl := simplify (polyl) : 
polyl := collect (polyl , h) : 

kd := alpha~3/(2-d*t)*coeff (polyl, h"n) : 
kd := normal (kd): 
kdn := kd: 

# 

# Extra manipulation for noncritical bundles with Omega=the Chern polynomial. 
# 

if s_diff_0(d) <> then 
kdn := kd/(s_dif f _0(d) ! ) : 
kd := x"(s_diff_0(d))*kd/(s_diff_0(d) !) : 
kd := normal (kd) 

fi: 

print ("kdn = " , kdn) : 
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print ("kd = ", kd) : 

# 

k := [ op(k), kd ] : 
print ("k = ", k) : 

# 

# Computation of $K_d$ via $Q_d$, using a formula in 

# Cox and Katz: Mirror symmetry and algebraic geometry 
# 

# 

# kda := alpha~2*coef f ( h*polyl, h~n )/d: 

# kda := normal (kda): 

# ka := [ op(ka) , kda ] : 

# print ("ka = ", ka) : 



# 

# Computation of instanton numbers from intersection numbers. 
# 

InstantonN := proc( inst_list, k_d ) 
local i, 1, nd, new_inst_list , s: 



1 := nops( inst_list ): 
nd := k_d: 



for i from 2 to 1+1 do 

if irem(l+l, i)=0 then 
s := iquo(l+l, i) : 
nd := nd - inst_list [s] /i~3 : 

f i: 

od: 



new_inst_list := [ op(inst_list) , nd ]: 

unassign( 'i', '1', 'nd', 's' ): 

RETURN(new_inst_list) : 

end: 



instanton_list := InstantonN( instanton_list , kdn) : 
print ("instanton_list = ", instanton_list) : 



od: 



unassign( 'd' , 'eqn_set', ' inhomogenll ' , ' inhomogenl2 ' , 'inhomogenl3' , 
' inhomogenl ' , ' inhomogen3 ' , ' inhomogen31 ' , ' inhomogen32 ' , 
'kl>, 'kla', 'kd', 'kda', 'kdn', '11', '12', 'nl', 'poly', 'polyl', 
'poly_set', 'poly_setl', 'poly_setll ' , 'poly_set21 ' , 'poly_set2', 
'poly_set3', 'poly_set4', ' qO_inverse ' , 'ql', 'ql_diff erence' , 
'qd', 'qdll', 'qdlll', 'qdll2', 'qdll3' , 'qdll4' , 'qdl2', 'qdl3', 
'qd21', 'qd22', 'qd23', 'qd31', 'qd311', 'qd312', 'qd313', 'qd32', 
'qdr', 'qdrl', 'qd_value', 'qd_valuel', 'qd_value2', 'qkl', 
'qkl_linear ' , 'rel', 'si', 's2', 's_diff, 'screen', 'substitution', 
' variable_list ' ): 



RETURN ( [q , k, instanton_list] ) : 
# RETURN ( [q , k, ka, instanton_list] ) : 
end: 



#■ 
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# SEC. 3 TEST ROUTINES and CASES OF STUDY. 

# EulerData(n, splitting_type , Omega_0, s_diff_0, dmax, opt_l, opt_2, check) 
# 

# REMARK. Machine used: Dell PC, Hard Drive 7 GB, 

# with Pentium II 400 MH, Cache 512 KB, and RAM 128 MB, upgrated to 384 MB on Dec. 26, 1999. 
# 

## 

## INSTRUCTION: For a bundle V of splitting type [[1_1, 1_2, ...],[k_l, k_2, ...]] over CP~n : 



## 

## (1) Omega: the Chern polynomial of the bundle V over CP~n. The indeterminate is denoted by x 
## and the hyperplane class of CP~n is denoted by h 

## (2) s_diff: the diffrence rank (U_d) -dim (M_{0,0}(CP~n, d)); see Sec. 2 of these notes for details. 
## it is a FUNCTION OF d . 

## E.g.\ 's_diff := d -> 3*d' is the function in usual notation ' s_dif f (d)=3d' . 

## (3) x : the indeterminate of the Chern polynomial. 

## x HAS TO BE SET TO ZERO IF s_diff IS A CONSTANT ZERO FUNCTION. 

## (Compare: Case 1 and Case 11.) 

## (4) opt_l : set to if only K_d and n_d are required; set to 1 if exact Q_d is required. 

## (5) opt_2 : if set to then 'LSolveSub' will be employed; if set to 1 then 'LSolveSub_A' will be 

## employed. Either is a subroutine to solve a system of linear equations. 

## (6) check : set to 1 if one wants to make sure that the system of linear equation is truely 

## consistent, thoughthe theory says it must be. Otherwise set to 0. 

## (7) dmax : the maximal degree one wants to compute. 

## (8) specialization: take Case 1 for example; the line 'lambda := i -> (i~2+7*i+l)*u' can be 

## replaced by e.g. 'lambda := i -> (i"3+17*i+2) *u' or even 'lambda := i-> i"3*u' 

## as long as making sure that at least visually that none of lambda(i)-lambda(j) for 

## $0\le i<j \le n$ is zero and that, when u=0, all lambda(i) must be zero. 

## Usually more complicated specialization has to be used for bigger dmax value. 

## (9) the first two argument of EulerData: n and splittingtype in the form [[...],[...]]. 

## Browse through the given examples below and the pattern should be clear. 

## 

## OUTPUT: To distinguish the two sets of numbers, the output contains built-in indications of 
## intersection numbers, "kd" , and instanton numbers, "instanton_list" . 

## 



## FOR NON-MAPLE-USER: 

## Note that the symbol '#' de-activates all the commands or words of the same line behind this symbol. 

## If you run the code as put here, it will compute only the case Case 1 since the command lines for 

## all other cases are de-activated by putting a '#' at the start. (Notice the difference with Mathematica) 

## To activate other cases, simply remove these '#' at the start of lines. 

## 

## Note, however, some of the lines are really meant to be titles or remarks and hence should be kept 
## de-activated; please compare Case 1 to keep these lines de-activated. 
## 

## CAUTION: While the numerals can be changed, all the variable names have to remain as given here 
## for the compatibility with the code. 

## 



## 

## Case 1 : 0(5) -> CP"4 

## 

restart : 

print ("Case 1: 0(5) -> CP'4."): 
x+5*h: 



s_diff 


= d 


-> 0: 








X 


= 0: 










opt_l 


= 0: 


###### 





or 


1 


opt_2 


= 0: 


###### 





or 


1 


check 


= 1: 


###### 





or 


1 


dmax 


= 3: 
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# specialization: 
# 

if opt_l=0 then 

lambda := i -> (i~2+7*i+l) *u 

f i: 

# 

qk := EulerData( 4, [[5], □], Omega, s_dif f , dmax, opt_l, opt_2, check ): 

print ( "Answer known n_d: 2875, 609250, 317206375, 242467530000, 229305888887625. 

print( "Case: QK(4, [[5], []]) = [", qk[2] , ", ", qk[3], "] . " ): 

# print (qk): ####### Re: This is de-activated for the output could be huge. 
unassign( 'check', 'dmax', 'lambda', 'Omega', 'opt_l', 'opt_2', 'qk', 

's_diff, 'x' ): 

## 



## 



Date: 


November 25, 1999. 


Specialization: 


lambda := i -> (i"4+3*i+l) *u 


Answer obtained K_d (d<=3) : 


2875, 4876875/8, 8564575000/27. 


K_d converted to n_d: 


2875, 609250, 317206375. 


Time consumed: 


52,034.1 sec. 



## 

## Case 2: 0(2)\oplus(4) -> CP"5 
## 

# restart: 

# print ("Case 2: 0(2) +0(4) -> 
x~2+6*h*x+8*h~2: 
d -> 0: 


###### or 
###### or 
###### or 



CP~5."): 



Omega 
s_diff 



# 
# 
# 
# 
# 
# 
# 

## 

## specialization: 
## 

# if opt_l = then 

# lambda := i -> (i"2+7*i+l)*u 

# fi: 



opt_l 
opt_2 
check 
dmax 



## 

# qk := EulerData( 5, [[2,4], [] ] , 

# print ( "Answer known n_d: 1280, 92288." ): 

# print( "Case: QK(5, [[2,4] , []] ) = [" , qk[2] 
## print (qk) : 

# unassign( 'check', 'dmax', 'lambda', 'Omega 
## 



s_diff, dmax, opt_l, opt_2, check ): 
, ", qk[3], "]." ): 
'opt_l', 'opt_2', 'qk', 's_diff 



## 



Date: 


November 24, 1999. 


Specialization: 


lambda := i-> (i~2+l)*u 


Answer obtained K_d (d<=2) : 


1280, 92448. 


K_d converted to n_d: 


1280, 92288. 


Consistency check: 


consistent . 


Time consumed: 


6,337.1 sec. 



## 

## Case 3: 0(3)\oplus 0(3) -> CP"5 
## 

# restart: 

# print ("Case 3: 0(3)\oplus 0(3) -> Cp-5."): 

# Omega := x~2+6*h*x+9*h~2 : 

# s_diff := d -> 0: 
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# X 


= 










# opt_l 


= 


###### 





or 


1 


# opt_2 


= 


###### 





or 


1 


# check 


= 1 


###### 





or 


1 


# dmax 


= 2 










## 













## specialization: 
## 

# if opt_l = then 

# lambda := i -> (i~2+l)*u 

# fi: 

## 

# qk := EulerData( 5, [[3,3], [] ] , Omega, s_dif f , dmax, opt_l, opt_2, check ): 

# print( "Answer known n_d: 1053, 52812." ): 

# print ( "Case: QK(5, [[3,3] , [] ] ) = [", qk[2] , ", ", qk[3], "] . " ): 
## print (qk) : 

# unassign( 'check', 'dmax', 'lambda', 'Omega', 'opt_l', 'opt_2', 'qk', 's_diff, 'x'): 



Date: 


November 24, 1999. 


Specialization: 


lambda := i -> (i"2+l)*u 


Answer obtained (d<=2) : 


1053, 423549/8. 


K_d converted to n_d: 


1053, 52812. 


Consistency check: 


consistent . 


Time consumed: 


5,693.4 sec. 



## 

## Case 4: 0(2)\oplus 0(2)\oplus 0(3) -> CP~6 
## 

# restart: 

# print ("Case 4: 0(2)\oplus 0(2) \oplus 0(3) -> CP"6."): 

# Omega := x~3+7*h*x"2+16*h~2*x+12*h~3 : 

# s_diff := d -> 0: 

# x := 0: 

# opt_l := 0: ###### or 1 

# opt_2 := 0: ###### or 1 

# check := 1: ###### or 1 

# dmax := 2: 

## 

## specialization: 
## 

# if opt_l = then 

# lambda := i -> (i~2+l)*u 

# fi: 

## 

# qk := EulerData( 6, [[2,2,3], [] ] , Omega, s_diff, dmax, opt_l, opt_2, check ): 

# print ( "Answer known n_d: 720, 22428." ): 

# print ( "Case: QK(6, [ [2, 2,3] , [] ] ) = [", qk[2], " , ", qk[3], "]." ): 
## print (qk) : 

# unassign( 'check', 'dmax', 'lambda', 'Omega', 'opt_l', 'opt_2', 'qk', 's_diff, 'x' ): 



Date: 


November 24, 1999. 


Specialization: 


lambda := i-> (i~2+l)*u 


Answer obtained K_d (d<=2) : 


720, 22518. 


K_d converted to n_d: 


720, 22428. 


Consistency check: 


consistent . 


Time consumed: 


50,205.2 sec. 



## 
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## Case 5: 0(2)\oplus 0(2) \oplus 0(2) \oplus 0(2) -> CP~7 
## 

# restart: 

# print ("Case 5: 0(2)\oplus 0(2)\oplus 0(2)\oplus 0(2) -> CP~7."): 



= x~4+8*h*x~3+24*h~2*x~2+32*h~3*x+16*h~4: 
= d -> 0: 
= 0: 

###### or 1 
###### or 1 
###### or 1 



= 

= 

= 1 

= 2 



# Omega 

# s_diff 

# x 

# opt_l 

# opt_2 

# check 

# dmax 

## 

## specialization: 
## 

# if opt_l = then 

# lambda := i -> (i+l)*u 

# fi: 

## 

# qk := EulerData( 7, [[2,2,2,2], []], Omega, s_diff, dmax, opt_l, opt_2, check ): 

# print( "Answer known n_d: 512, 9728." ): 

# print ( "Case: QK(7, [ [2, 2, 2, 2] , [] ] ) = [", qk[2], ", ", qk[3], "]." ): 
## print (qk) : 

# unassign( 'check', 'dmax', 'lambda', 'Omega', 'opt_l', 'opt_2', 'qk', 's_diff, 'x' ): 
## 

## Date: January 8, 2000 

## Specialization: lambda := i -> (i"2+13*i+l) *u 

## Answer obtained K_d (d=l) : 512, . 

## K_d converted to n_d : 512. 

## Consistency check: . 

## Time consumed: 91.3 sec 

## Bytes used: 22.7 MB 

## 

## Remark. The computaion of K_2 exceeds the capacity of Maple V, Release 5.1. 
## 



## 

## Case 6: T_{\ast}CP"2 
## 

# restart: 

# print ("Case 6: T_{\ast}CP~2 . ") : 



= x~2+3*h*x+3*h~2: 
= d -> 3: 



= 

= 1 

= 1 

= 3 



###### or 1 
###### or 1 
###### or 1 



# Omega 

# s.diff 

# opt_l 

# opt_2 

# check 

# dmax 

## 

## specialization: 
## 

# if opt_l = then 

# lambda := i -> (i"3+17*i+l) *u 

# fi: 



# qk := EulerData( 2, [[1,2], [] ] , Omega, s_dif f , dmax, opt_l, opt_2, check ): 

# print ( "Answer known: Unknown."): 

# print ( "Case: QK(2, [[1,2] , []] ) = [" , qk[2] , ", ", qk[3], "] . " ): 
## print ("qk = ", qk) : 

# unassign( 'check', 'dmax', 'lambda', 'Omega', 'opt_l', 'opt_2', 'qk', 

# 's_diff, 'x' ): 

## 

## Date: December 8, 1999. 
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## Specialization: lambda := i -> (i~2+3*i+l)*u 

## Answer obtained K_d (d=l) : 10*x~3 

## Consistency check: ??? 

## Time consumed: ??? sec. 

## Bytes consumed: ??? MB 

## Remark: For the computation of d=2, the complexity involved seems to surpass the capacity of Maple V. 

## 

## Date: January 7, 2000. 

## Specialization: lambda := i -> (i~3+17*i+l) *u 

## Time consumed: ??? 

## Bytes consumed: ??? > 150 MB 

## 

## Remarks: So far, the job has been thrown out constantly at d=2. After pinning down exact what line 



## that cause an error message to appear and the job thrown out, it is found that, in the 

## computation for d=2, a coeffiecent (which is a rational function of variable 'x' and 'u') 

## that takes nearly 900 lines for the display in the Maple window appears. This seems to beat 

## the current capacity of Maple. Finding a way to go around this will be one of the major 

## improvements of the current code. 

## 



## 

## Case 7: 0(-l)\oplus 0(-l) -> CP~1 
## 

# restart: 

# printC'Case 7: 0(-l)\oplus 0(-l) -> CP~1"): 



# Omega := l/(x-h)~2: 

# s_diff := d -> 0: 

# x := 0: 

# opt_l := 0: ###### or 1 

# opt_2 := 0: ###### or 1 

# check := 1: ###### or 1 

# dmax := 5: 

## 



## specialization: 
## 

# if opt_l = then 

# lambda := i -> (i+l)*u 

# fi: 

## 

# qk := EulerData( 1, [ [] , [-1, -1]], Omega, s_diff, dmax, opt_l, opt_2, check ): 

# print( "Answer known: K_d = l/(d"3)." ): 

# print ( "Case: QK(1, [ [] , [-1 , -1] ] ) = [", qk[2], ", ", qk[3], "]." ): 
## print (qk) : 

# unassign( 'check', 'dmax', 'lambda', 'Omega', 'opt_l', 'opt_2', 'qk', 

# 's_diff, 'x' ): 

## 

## $K_d$ tested correct for d<=5. 

## 

## 

## Case 8: 0(-3) -> CP"2 

## 

# restart: 

# printC'Case 8: 0(-3) -> CP"2") : 



# Omega := l/(x-3*h) : 

# s_diff := d -> 0: 

# x := 0: 

# opt_l := 0: ###### or 1 

# opt_2 := 0: ###### or 1 

# check := 1: ###### or 1 
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# dmax := 5: 

## 

## specialization: 
## 

# if opt_l = then 

# lambda := i -> (i~3+i~2+3*i+l)*u 

# fi: 

## 

# qk := EulerData( 2, [ [] , [-3]], Omega, s_diff, dmax, opt_l, opt_2, check ): 

# print ( "Answer known K_d: 3, -45/8, 244/9, -12333/64, 211878/125." ): 

# print ( "Case: QK(2, [[],[-3]]) = [", qk[2], ", ", qk[3], "] . " ): 
## print (qk) : 

# unassign( 'check', 'dmax', 'lambda', 'Omega', 'opt_l', 'opt_2', 'qk', 's_diff, 'x' ): 
## 

## Date: November 26, 1999. 

## Specialization: lambda := i -> (i~3+i~2+3*i+l) *u 

## Answer obtained K_d (d<=5) : same as the known ones. 

## Consistency check: consistent. 

## Time consumed: 48,120 sec. 

## 



## 

## Case 9: 0(2)\oplus 0(-2) -> CP"3 
## 

# restart: 

# printC'Case 9: 0(2)\oplus 0(-2) -> CP~3."): 



= (x+2*h)/(x-2*h) : 
= d -> 0: 
= 0: 

###### or 1 
###### or 1 
###### or 1 



= 

= 

= 1 

= 3 



# Omega 

# s.diff 

# x 

# opt_l 

# opt_2 

# check 

# dmax 

## 

## specialization: 
## 

# if opt_l = then 

# lambda := i -> (i~2+7*i+l)*u 

# fi: 

## 

# qk := EulerData( 3, [[2], [-2]], Omega, s_diff, dmax, opt_l, check ): 

# print ( "Answer known K_d: -4, -9/2, -328/27." ): 

# print ( "Case: QK(3, [[2], [-2]]) = [», qk[2] , ", ", qk[3], "] . " ): 
## print (qk) : 

# unassign( 'check', 'dmax', 'lambda', 'Omega', 'opt_l', 'opt_2', 'qk', 's_diff, 'x' ): 
## 

## $K_d$ tested correct for d<=3. 

## 



## 

## Case 10: 0(2)\oplus 0(2)\oplus 0(-l) -> CP~4 
## 

# restart: 

# printC'Case 10: 0(2)\oplus 0(2)\oplus 0(-l) -> CP"4."): 



# Omega 

# s_diff 

# x 

# opt_l 

# opt_2 

# check 



= (x+2*h)"2/(x-h) : 
= d -> 0: 
= 0: 

= 0: ###### or 1 
= 0: ###### or 1 
= 1: ###### or 1 
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# dmax := 2: 

## 

## specialization: 
## 

# if opt_l = then 

# lambda := i -> (i~2+l)*u 

# fi: 

## 

# qk := EulerData( 4, [[2,2], [-1]], Omega, s_diff, dmax, opt_l, check ): 

# print( "Answer known K_d: 16, -18, 1312/27." ): 

# print ( "Case: QK(4, [ [2, 2] , [-1] ] ) = [", qk[2], ", ", qk[3], qk[4], "] . " ): 
## print (qk) : 

# unassign( 'check', 'dmax', 'lambda', 'Omega', 'opt_l', 'opt_2', 'qk', 's_diff, 'x' ): 
## 

## $K_d$ tested correct for d<=2 

## 



## 

## Case 11: 0(2) -> CP"1 
## 

# restart: 

printC'Case 11: 0(2) 
x+2*h: 



-> CP"1.") : 



s_diff 

opt_l 

opt_2 

check 

dmax 



-> 3: 



###### or 
###### or 
###### or 



# 
# 
# 
# 
# 
# 
# 

## 

## specialization: 
## 

# if opt_l = then 

# lambda := i -> (i~2+l)*u 

# fi: 

## 

# qk := EulerData( 1, [[2], []], 

# print ( "Case: QK(1, [[2], []]) 
## print (qk) : 

# unassign( 'check', 'dmax', 
## 

## Date: 

## Specialization: 
## Answer obtained K_d (d<=5) 
## K_d converted to n_d 
## Consistency check: 
## Time consumed: 
## Bytes used: 

## 



[", 



s_diff, dmax, opt_l, opt_2, check ): 
qk[2], ", ", qk[3], "] . » ): 



lambda', 'Omega', 'opt_l', 'opt_2', 'qk', 's_diff 



December 8, 1999. 
lambda := i -> (i~2+l)*u 

x"3, (l/8)*x~3, (l/27)*x~3, (l/64)*x"3, (l/125)*x"3 
x"3, 0, 0, 0, 0. 
consistent . 
134.5 sec. 
3.25 MB 



## Case 12: 0(3) -> CP"2 

## 

# restart: 

# printC'Case 12: 0(3) -> Cp-2."): 



# 


Omega 


= x+3*h 










# 


s_diff 


= d -> 


2: 








# 


opt_l 


= 0: 


###### 





or 


1 


# 


opt_2 


= 0: 


###### 





or 


1 


# 


check 


= 1: 


###### 





or 


1 


# 


dmax 


= 6: 
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## 

## specialization: 
## 

# if opt_l = then 

# lambda := i -> (i~2+7*i+l) *u 

# fi: 



## 

# qk := EulerData( 2, [[3], 

# print ( "Case: QK(2, [[3], 
## print (qk) : 

# unassign( 'check', 'dmax' 
## 



[]], Omega, s_diff, 
[]]) = [", qk[2], " 



dmax, opt_l, opt_2, check ): 
, ", qk[3], »].» ): 



'lambda' 



' Omega' , ' opt_l ' 



'opt_2' 



'qk' , 



's_diff ' 



???, (Recording date: December 23, 
lambda := i -> (i~2+7*i+l) *u 



1999) . 



## Date: 

## Specialization: 
## Answer obtained K_d (d<=6) : 

## 21*x"2, (189/8)*x~2, (169/9)*x~2, (1533/64) *x~2 , (2646/125)*x~2, 

## K_d converted to n_d: 21*x"2, 21*x~2, 18*x~2, 21*x~2, 21*x~2, 

## Consistency check: consistent. 

## Time consumed: 1,047,896.9 sec. 

## Bytes used: 116 MB 

## 



(169/8) *x"2 
18*x~2, 



restart : 

print ("Case 13: 



## 

## Case 13: 0(4) -> CP"3 
## 

# 
# 
# 
# 
# 
# 
# 
# 

##- 



0(4) -> CP"3."): 



Omega 

s_diff 

opt_l 

opt_2 

check 

dmax 



x+4*h: 
1: 



1 

5 



###### or 
###### o or 
###### o 



1 
1 

or 1 



## specialization: 
## 

# if opt_l = then 

# lambda := i -> (i"2+7*i+l)*u 

# fi: 

## 



# qk := EulerData( 3, [[4], []], Omega, s_diff, dmax, opt 

# print ( "Case: QK(3, [[4], []]) = [", qk[2], ", ", qk[3] 
## print (qk) : 

# unassign( 'check', 'dmax', 'lambda', 'Omega', 'opt_l', 'opt_2' 
## 



1 , opt_2, check ) : 
, "]." ): 



'qk' 



's_diff ' 



## Date: 

## Specialization: 

## Answer obtained K_d (d<=4) 

## K_d converted to n_d 

## Consistency check: 

## Time consumed: 

## Bytes used: 

## 



December 8, 1999. 

lambda := i -> (i~2+7*i+l) *u 

320*x, 5056*x, (3893504/27) *x, 5490624*x, 

320*x, 5016*x, 144192*x, 5489992*x. 

consistent . 

55,290.4 sec 

30.2 M 



## 

## Case 14: 0(6) -> CP"4 

## 

# restart: 

# printC'Case 14: 0(6) -> CP"4") : 

# Omega := x+6*h: 
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d -> 

0: 

1: 

1: 

10: 



# s.diff 

# opt_l 

# opt_2 

# check 

# dmax 

## 

## specialization: 
## 

# if opt_l =0 then 

# lambda := i -> 

# fi: 

## 

# qk 



d: 
###### 
###### 
###### 



(i~2+7*i+l)*u 



[[6], 
[[6], 

' dmax ' , 



:= EulerData( 4, 

# print ( "Case: QK(4, 
## print (qk) : 

# unassign( 'check', 
## 

## Date: 

## Specialization: 

## Answer obtained K_d (d<=3) 

## K_d converted to n_d 

## Consistency check: 

## Time consumed: 

## Bytes used: 

## 



□ ]. 
[]]) 



s_diff, dmax, opt_l, opt_2, 
qk[2], », », qk[3], »] . " ): 



check ) : 



'lambda', 'Omega', 'opt_l', 'opt_2', 'qk', 's_diff', 'x' ): 



December 28, 1999. 

lambda := i -> (i~2+7*i+l) *u 

50400*x, (752729895/4) *x~2, (433244745198080/243) *x~3 , — 
50400, 752704695/4, 433244744744480/243, 
consistent . 
49,842.4 sec 
30.5 M 



15: 0(7) -> CP"4 



restart : 

print ("Case 15: 0(7) -> CP~4."): 
x+7*h: 
d -> 2*d: 
0: 



## 

## Case 

## 

# 
# 
# 
# 
# 
# 
# 
# 

## 

## specialization: 
## 

# if opt_l = then 
lambda : 



Omega 
s.diff 
opt_l 
opt_2 
check 
dmax 



0: 
1: 
10: 



###### 
###### 
###### 



> (i"2+ll*i+l)*u 



[[7], 
[[7], 

' dmax ' 



# 

# fi: 

## 

# qk := EulerData( 4 

# print ( "Case: QK(4 
## print (qk) : 

# unassign( 'check', 
## 

## Date: 

## Specialization: 

## Answer obtained K_d (d<=3) 

## K_d converted to n_d : 

## Consistency check: 

## Time consumed: 

## Bytes used: 

## 



[]] 

[]]) = [" 
'lambda' 



s_diff , 
qk[2], 1 



dmax, opt_l, opt_2, 
, ", qk[3] , »].» ): 



check ) : 



'opt_l', 'opt_2', 'qk', 's.diff 



): 



December 29, 1999. 

lambda := i -> (i~2+ll*i+l) *u 

451570*x~2 , (403985396325/32) *x"4 , (15755269694706695755/17496) *x~6 , 
451570, 403983590045/32, 15755269694414078395/17496, — 
consistent . 
72,086.1 sec. 
45.2 MB. 



## 

## Case 16: 0(8) -> CP"4 
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## 

# restart: 

print ("Case 16: 0(8) 
x+8*h: 



-> CP-4.") : 



s_diff 
opt_l 
opt_2 
check 
dm ax 



d -> 

0: 

0: 

1: 

10: 



3*d: 
###### or 
###### or 
###### or 



(i~2+ll*i+l)*u 



# 
# 
# 
# 
# 
# 
# 

## 

## specialization: 
## 

# if opt_l = then 

# lambda := i -> 

# fi: 

## 

# qk := EulerData( 4, [[8], []], Omega, s_dif f , 

# print( "Case: QK(4, [[8], []]) = [", qk[2], " 
## print (qk) : 

# unassign( 'check', 'dmax', 'lambda', 'Omega', 
## 

## Date: December 30, 1999. 

## Specialization: lambda := i -> (i~2+ll*i+l) *u 

## Answer obtained K_d (d<=3) : 

## 2773820*x~3, (3178734062035/8) *x"6 , (46028387589557254161275/314928) *x"9 , — . 

## K_d converted to n_d : 2773820, 3178731288215/8, 46028387589524900324795/314928, — 

## Consistency check: consistent. 

## Time consumed: 80,553.2 sec. 

## Bytes used: 49.8 MB. 

## 



dmax, opt_l, opt_2, check ): 
, ", qk[3] , "]." ): 

'opt_l', 'opt_2', 'qk', 's_diff', 'x' ): 



## 

## Case 17: 0(9) -> CP"4 

## 

# restart: 

# print ("Case 17: 0(9) -> CP"4."): 

# Omega := x+9*h: 

# s_diff := d -> 4*d: 

# opt_l := 0: ###### or 1 

# opt_2 := 0: ###### or 1 

# check := 1: ###### or 1 

# dmax := 10: 

## 

## specialization: 
## 

# if opt_l = then 

# lambda := i -> (i"2+ll*i+l) *u 

# fi: 

## 

# qk := EulerData( 4, [[9], []], Omega, s_diff, dmax, opt_l, opt_2, check ): 

# print( "Case: QK(4, [[9], []]) = [", qk[2], ", ", qk[3], "] . » ): 
## print (qk) : 

# unassign( 'check', 'dmax', 'lambda', 'Omega', 'opt_l', 'opt_2', 'qk', 's_diff, 'x' ): 
## 

## Date: December 31, 1999. 

## Specialization: lambda := i -> (i"2+ll*i+l) *u 

## Answer obtained K_d (d<=3) : 

## 13198850*x"4, (243281907041715/32) *x"8 , (197802281929974511821535/17496)*x-12, — . 

## K_d converted to n_d : 13198850, 243281854246315/32, 197802281929965958966735/17496, 

## Consistency check: consistent. 

## Time consumed: 93,665.7 sec. 
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## Bytes used: 
## 



53.2 MB. 



## 

## Case 18: 0(10) -> CP~4 
## 

# restart: 

# print ("Case 18: 0(10) -> CP'4."): 

# Omega := x+10*h: 

# s_diff := d -> 5*d: 

# opt_l := 0: ###### or 1 

# opt_2 := 0: ###### or 1 

# check := 1: ###### or 1 

# dmax := 10: 

## 

## specialization: 
## 

# if opt_l = then 

# lambda := i -> (i~2+ll*i+l) *u 

# fi: 

## 

# qk := EulerData( 4, [[10], [] ] , Omega, s_dif f , 

# print ( "Case: QK(4, [[10], [] ] ) = [" , qk[2], " 
## print (qk) : 

# unassign( 'check', 'dmax', 'lambda', 'Omej 

## 

## Date: January 2, 2000. 

## Specialization: lambda := i -> (i~2+ll*i+l) *u 

## Answer obtained K_d (d<=3) : 

## 52040450*x"5, (25908993204089625/256) *x"10 , (71418501571607082433686025/139968) *x~15 , — . 
## K_d converted to n_d : 52040450, 25908991538795225/256, 71418501571606812655993225/139968, 

## Consistency check: consistent. 
## Time consumed: 80,874.0 sec. 

## Bytes used: 49.6 MB. 

# 



dmax, opt_l, opt_2, check ): 
, ", qk[3], "]." ): 



'opt_l', 'opt_2' 



'qk' , 



's_diff ' 



): 



End of the code . 
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